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Abstract. The commonly used classical equipartition or minimum-energy estimate of total magnetic fields strengths from 
radio synchrotron intensities is of limited practical use because it is based on the hardly known ratio K, of the total energies 
of cosmic ray protons and electrons and also has inherent problems. We present a revised formula, using the number density 
ratio K for which we give estimates. For particle acceleration in strong shocks K is about 40 and increases with decreasing 
shock strength. Our revised estimate for the field strength gives larger values than the classical estimate for flat radio spectra 
with spectral indices of about 0.5-0.6, but smaller values for steep spectra and total fields stronger than about 10 /iG. In 
very young supernova remnants, for example, the classical estimate may be too large by up to lOx. On the other hand, 
if energy losses of cosmic ray electrons are important, K increases with particle energy and the equipartition field may be 
underestimated significantly. Our revised larger equipartition estimates in galaxy clusters and radio lobes are consistent with 
independent estimates from Faraday rotation measures, while estimates from the ratio between radio synchrotron and X-ray 
inverse Compton intensities generally give much weaker fields. This may be explained e.g. by a concentration of the field in 
filaments. Our revised field strengths may also lead to major revisions of electron lifetimes in jets and radio lobes estimated 
from the synchrotron break frequency in the radio spectrum. 

Key words: ISM: magnetic fields - ISM: supernova remnants - galaxies: active - galaxies: clusters: general - galaxies: 
magnetic fields - radio continuum: general 
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1. Introduction 

Magnetic fields play an important role in supernova remnants 
(Fiirst & Reich 2004\ in the interstellar medium of galaxies 
(Beck et al. Beck I20()4t . in the intergalactic medium 

of clusters (Carilli & Tavlor l2002> . and in the jets and lobes 
of active galaxies (Rees [TW51 Blandford 200X1. The relative 
importance among the other competing forces can be esti- 
mated by comparing the corresponding energy densities or 
pressures. For most extragalactic objects measurements of 
the magnetic field strength are based on radio synchrotron 
emission and Faraday rotation measures (RM) of the polar- 
ized radio emission (see Heiles [T976l and Verschuur l 1 979l for 
reviews of the observational methods). 

Total synchrotron emission traces the total field in the sky 
plane, polarized synchrotron emission the regular field com- 
ponent. The regular field can be coherent (i.e. preserving its 
direction within the region of observation) or incoherent (i.e. 
frequently reversing its direction). 
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Determination of the field strength from the synchrotron 
intensity needs information about the number density of the 
cosmic ray electrons, e.g. via their 7-ray bremsstrahlung 
emission or X-ray emission by inverse Compton scattering. 
If such data are lacking, an assumption about the relation 
between cosmic ray electrons and magnetic fields has to be 
made. The most commonly used approaches are: 

- Minimum total energy density (etot — ecR + eb = min) 

- Equipartition between the total energy densities of cosmic 
rays and that of the magnetic field (ee = ecR) 

- Pressure equilibrium (eB = lecR) 

The minimum-energy and equipartition estimates give very 
similar results and are often used synonymously. The idea is 
that cosmic ray particles and magnetic fields are strongly cou- 
pled and exchange energy until equilibrium is reached. Devi- 
ations from equilibrium occur if, for example, the injection 
of particles or the generation of magnetic fields is restricted 
to a few sources or to a limited period. The mean distance 
between sources and their mean lifetime define the smallest 
scales in space and time over which equipartition can be as- 
sumed to hold. 



2 THE CLASSICAL MINIMUM-ENERGY FORMULA 



The minimum-energy assumption was first proposed by 
Burbidge (1956 ) and applied to the optical synchrotron emis- 
sion of the jet in M87. Since then the vahdity of this method 
has been discussed in the literature. Duric (I1990> argued that 
any major deviation from equipartition would be in conflict 
with radio observations of spiral galaxies. 

The azimuthally averaged equipartition strength of the 
field in the Milky Way and its radial variation (Berkhui- 
jsen, in Beck 1200 H agree well with the model of Strong 
et al. (IJOOO their Fig. 6) based on radio continuum and 7- 
ray surveys and cosmic ray data near the sun. On the other 
hand, Chi & Wolfendale (1993 1 claimed significant devia- 
tions from equipartition conditions in the Magellanic Clouds. 
Pohl I ll993b> replied that the standard proton-to-electron ra- 
tio used in the equipartition formula may be smaller in the 
Magellanic Clouds compared to the Milky Way. 

Equipartition estimates of field strengths were determined 
in many spiral galaxies (Beck 2000 1. The mean strength of 
the total field ranges from a few in galaxies with a low 
star-formation rate to ~ 30 in grand-design galaxies like 
M 5 1 . The mean total field strength of a sample of Shapley- 
Ames galaxies is 9 piG (Niklas 1997 see Sect. l4.lt . The re- 
lation between the total field strength and the star formation 
rate is deduced from the radio-infrared correlation (Helou & 
Bicay|l993 Niklas & Beck 1997 1. The total field is strongest 
in the spiral arms, but the strength of the large-scale regu- 
lar field in most galaxies is highest in interarm regions (Beck 
I2OOOII2OOII . 

In the radio lobes of strong FRII radio galaxies the 
equipartition estimates of the field strength are 20-100 /iG, 
in hot spots 100-600 /iG, and 0.5-1 /iG in the intracluster 
medium of galaxy clusters and in relics (see Sect. 14. 3> . In 
these estimates relativistic protons were assumed to be ab- 
sent ("absolute minimum energy") or contribute similarly to 
the total energy density as the electrons. The ratio of radio 
synchrotron to inverse Compton X-ray intensities can be used 
as another estimate of the field strength (e.g. Carilli & Taylor 
I2002> . In most radio lobes the two estimates are similar, but 
there are also significant deviations where "out of equiparti- 
tion" conditions have been claimed (see Sect. l4.4> . 

Faraday rotation measures RM are sensitive to the co- 
herent regular field component along the line of sight and 
to the density of thermal electrons which can be derived 
from thermal emission or pulsar dispersion measures DM. 
The ratio RM/DM is a widely used measure of coherent 
field strengths in the Milky Way. The derived values are 
lower than the equipartition estimates from polarized inten- 
sity (Beck et al. 2003). In galaxy clusters, observations of 
fluctuations in RM are used to estimate total field strengths. In 
the Coma cluster such data indicate fields one order of mag- 
nitude stronger than the equipartition value (Giovannini et al. 
[T993I Fer etti et al. lT995l Fusco-Femiano et al. lT999l Govoni 
& Feretti "2004"! . Strong total fields have also been claimed 
from observations of rotation measures in many other clus- 
ters (CarilU & Tavlor l2002l . 

All methods to estimate field strengths are subject to bias. 
Firstly, a positive correlation between field strength and cos- 
mic ray electron density on small scales, which are unre- 



solved by the telescope beam or occur along the line of sight, 
leads to an overestimate of the equipartition strengths of the 
total and regular field components (Beck et al. 2003 1. Fur- 
thermore, synchrotron intensity is biased towards regions of 
strong fields, so that the equipartition estimates are higher 
than the average field strength (Sect. I4.4> . Field strengths 
based on Faraday rotation measures may also be overesti- 
mated if small-scale fluctuations in magnetic field and in ther- 
mal gas density are positively correlated. Finally, Newman et 
al. ('2002'l and Murgia et al. (2004 1 pointed out that field esti- 
mates based on Faraday rotation measures are likely to be too 
high by a factor of a few if there is a spectrum of turbulence 
scales. 

In this paper we show that there are inherent problems 
with the classical equipartition / minimum-energy formula, 
especially when computing its parameters from observable 
quantities. We present a revised formula which is based on 
integration of the cosmic ray energy spectrum rather than the 
radio frequency spectrum and discuss the limitations of its 
application. 

2. The classical minimum-energy formula 

The classical textbook minimum-energy formula is based on 
a value for the total cosmic ray energy density which is deter- 
mined by integrating the radio spectrum from i^min to I'max, 
usually from 10 MHz to 10 GHz, which is the typical range 
accessible to observations. Particles and magnetic fields are 
assumed to fill the source homogeneously, and the field is 
assumed to be completely tangled (no preferred direction). 
Then the field minimum-energy strength -Bmin is quoted as 
(e. g. Pach olcyk lT97()l Longair [T994l p. 292; Rohlfs & Wil- 
son [T996l : 

Selass= (6^G(^+1)L,/F)'/' (1) 

where G is a product of several functions varying with h'rnm, 
Vniax and synchrotron spectral index a (see Eq. (IA20> for 
details). K, is the ratio of the total energy of cosmic ray nuclei 
to that of the synchrotron emitting electrons + positrons. is 
the synchrotron luminosity at frequency ly, V is the source's 
volume, and L^/V is the synchrotron emissivity. 

The resulting magnetic energy density ee is 3/4 of the 
cosmic ray energy density ecR, so that -Bciass is 8% smaller 
than the equipartition field strength -Bcq,ciass: 

Sclass/-Bcq, class = (3/4)^^'' 

The problems with the classical formula are the follow- 
ing: 

(1) The radio emissivity L^/V is the average over the 
source's volume. Since radio intensity, in the case of equipar- 
tition between magnetic fields and cosmic rays, varies with 
about the fourth power of the total field strength, it emerges 
mainly from the regions with the highest field strengths, so 
that Smin is larger than the volume-averaged field strength if 
the local field strength B is inhomogeneous within the source. 
This problem can be reduced by replacing L^/V by ly/l. 
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where Ii^ is the local synchrotron intensity (surface bright- 
ness) and / is the pathlength through the emitting medium 
(see Eq. (lA20l l 

(2) A fixed integration interval Vaiin to Vnia^ is used in 
Eq. Q. The critical frequency of synchrotron emission of an 
electron with energy £' in a magnetic field of strength B with 
a component B± perpendicular to the velocity vector is (Lang 
11^9 p. 29): 

= ciE^BjL = 16.08 MHz (i;/GeV)2 (B_l/mG) (2) 

where ci — 3e/(47rmo^c^). e is the elementary charge, nic 
the electron mass, and c the velocity of light. 

Note that the synchrotron emission of a single electron is 
maximum at the frequency Vniax — 0.29i'crit (Longair l 19941 
p. 246), but for a continuous energy spectrum of electrons 
the maximum contribution at a given frequency is emitted by 
electrons which are a factor of almost two lower in energy 
(Webber et al. ll980t . As a result, Eq. (|2ji can be used for the 
maximum synchrotron emission from a spectrum of electrons 
around E. 

The standard integration interval of 10 MHz- 10 GHz cor- 
responds to an interval [Ei - i?2] in the electron energy 
spectrum of 800 MeV-25 GeV in a 1 /xG field, 250 MeV- 
8 GeV in a 10 ^G field, or to 80 MeV-2.5 GeV in a 100 iiG 
field. The integrated cosmic ray energy is proportional 
to [E\~^°' — i?2~^"] where a is the synchrotron spectral in- 
dex (see Eq. ( IA5» . and Ei and E2 are fixed integration lim- 
its. Replacing Ei and E2 by t'min and i^max in the classical 
minimum-energy formula via Eq. (|2j introduces an additional 
term depending on the magnetic field strength. As a conse- 
quence, the total energy etot depends on a constant power 
of B (see Longair [T994I p. 292) and the wrong derivative 
dctot/dB leads to the constant exponent 2/7 in Eq. Q. The 
classical minimum-energy formula is formally incorrect. 

(3) K, in Eq. Q is the ratio of the total energy density 
of cosmic ray nuclei to that of the electrons (H-positrons). 
Knowledge of K, would require measurements of the spectra 
of the main cosmic ray species over the whole energy range, 
especially at low particle energies which, for steep spectra, 
contribute mostly to the total energy. Instead, the total en- 
ergy of cosmic ray electrons is approximated in the classi- 
cal formula by integrating the radio spectrum between fixed 
frequency limits, and the energy of the cosmic ray nuclei is 
assumed to scale with JC. 

This classical procedure is subject to major uncertain- 
ties. Firstly, the observable radio spectrum only traces the 
spectrum of cosmic ray electrons over a small energy range 
(Eq. Secondly, the ratio K, of total energies may dif- 
fer from that in the observable energy range. What would 
be needed in the classical formula is the energy ratio K,' in 
the limited energy range traced by the observed synchrotron 
emission. In case of energy losses (Sect. lT4l i JC' may deviate 
strongly from JC. As the other input numbers of the classical 
formula are generally known with sufficient accuracy, the un- 
certainty in JC is the reason why the equipartition / minimum 
energy field strengths are regarded as crude estimates and the 
field estimates are often given with a scaUng of {JC + l)-^/''. 



We propose to use, instead of the energy ratio JC, the ra- 
tio of number densities K of cosmic ray protons and elec- 
trons per particle energy interval within the energy range 
traced by the observed synchrotron emission. Measurements 
of the local Galactic cosmic rays near the sun (see Appendix) 
yield Kq ~ 100 at a few GeV, which is the energy range 
relevant for radio synchrotron emission. This value is con- 
sistent with the predictions from Fermi shock acceleration 
and hadronic interaction models (see Table [Q. At lower and 
higher energies, however, K may vary with particle energy 
(see Sect.lsH. 

The observed energy spectrum of cosmic rays is the re- 
sult of balance between injection, propagation and energy 
losses. Interactions with matter and radiation are different 
for protons and electrons, so that the shape of their energy 
spectra are generally different (Pohl ll993al Lieu et al. 1 19991 
Schlickeiser '2002'!. At low energies (typically below a few 
100 MeV) the dominating loss of cosmic ray protons and 
electrons is ionization of the neutral gas and/or Coulomb in- 
teraction with the ionized gas. At energies beyond 1 GeV 
the dominating loss of protons (and other nucleons) are in- 
elastic collisions with atoms and molecules, producing pions 
and secondary electrons. The spectral index of the stationary 
energy spectrum is not changed. The dominating loss mech- 
anism for electrons is nonthermal bremsstrahlung producing 
low-energy 7-rays (Schlickeiser 2002 p. 100). At even higher 
energies (beyond a few GeV) the electrons suffer from syn- 
chrotron and inverse Compton losses (see Sect. I3.4> . Fur- 
thermore, the spectra of all cosmic ray species in galaxies 
may be steepened if particle diffusion (escape) is energy- 
dependent. As the result of energy loss processes, the cos- 
mic ray electron spectrum is not proportional to the proton 
spectrum, so that K varies with energy. Only in a relativis- 
tic electron/positron plasma as discussed for jets and lobes of 
radio galaxies, where cosmic ray protons are absent, JC = 
and K = are valid at all energies (see Sects. 13. 3l and l4.3> . 

In this paper we present an easily applicable formula 
with two input parameters from observations, synchrotron in- 
tensity and spectral index. We discuss the energy/frequency 
range where the revised formula can be applied because a re- 
liable and constant value for the proton-electron number den- 
sity ratio Kq can be adopted. As a result, a more accurate 
estimate of the equipartition field strength is possible than in 
the classical approach. 

Pfrommer & EnBlin (2004b') give two formulae, for the 
classical minimum-energy criterion and for the hadronic in- 
teraction model applied to two galaxy clusters, taking into ac- 
count (2) and (3) discussed above. Their formula (6.2) for the 
classical case includes luminosity, cluster volume, the proton- 
to-electron energy density ratio and the lower cut-off energy 
of the electron spectrum. 

3. The revised formula 

The equipartition and minimum-energy procedures need the 
total energy of cosmic rays, to be derived from the observed 
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3.2 Revised minimum energy formula for the total field 



radio synchrotron spectrum. An accurate treatment needs to 
account for all energy loss processes which affect the energy 
spectrum of nucleons and electrons, as discussed in detail 
by Pohl (I1993a> . This method, however, requires additional 
knowledge of the distributions of gas density and radiation 
field and is applicable only to a few well-studied objects. 

In this paper we derive a revised formula for the simple 
case that the number density ratio Kq of protons and elec- 
trons is constant which is valid in a limited range of particle 
energies. We further assume that the cosmic rays are accel- 
erated by electromagnetic processes which generate power 
laws in momentum, and that the same total numbers of pro- 
tons and electrons are accelerated. The total energy is domi- 
nated by the protons. As the proton energy spectrum flattens 
below the fixed proton rest energy Ep, the total cosmic ray 
energy can be computed easily. The details are presented in 
the Appendix. 

3.1. Revised equipartition formula for the total field 

From the total cosmic ray energy density ecR as a function 
of field strength B (Eq. jA16L see Appendix) and assuming 

ecR = cb = B'^J 8n we get: 



Beq = {4^(2a+l)(Ko + l)/. ^-2" W2cir 



/ [i2a-l)c2{a)lc4{i)]} 



1/(q+3) 



(3) 



is the synchrotron intensity at frequency ly and a the syn- 
chrotron spectral index. Kq is the constant ratio of the num- 
ber densities of protons and electrons in the energy range 
where their spectra are proportional (see Sect. l3.4> . Ep is the 
proton rest energy, ci, C2 and C4 are constants. C2 and C4 de- 
pend on a and the magnetic field's incUnation, respectively 
(see the Appendix for details). 

Ii, and a can be determined from observations, while the 
proton-to-electron ratio Kq and pathlength I have to be as- 
sumed. If the synchrotron sources have a volume filling fac- 
tor /, I has to be replaced by / x / in order to obtain B^q 
within the source. B^q depends only weakly on the source's 
distance via the dependence on I. 

We have restricted the discussion in this paper to nearby 
sources. For redshifts z > 0, correction terms are required 
which are given e.g. in Govoni & Feretti (2004). 

Eq. Q yields field strengths which are larger by 7% for 
a = 0.6, larger by 2% for a — 0.75, and smaller by 8% 
for a — 1.0 compared to the results obtained with the ear- 
lier version of the revised equipartition formula (e.g. Krause 
et al.[T984, Beck 1991, Niklas 1997; Thierbachet al. 2003). 
The reason is that a simplified version of Eq. ( IA12t was used 
previously, with integration only from Ei = 300 MeV to 
E2 ^ oo- However, the differences to Eq. Q are small, 
smaller than the typical errors caused by uncertainties in the 
input values, so that the values published by the Bonn group 
previously are still valid. 

In case of equilibrium between magnetic and cosmic ray 
pressures, the field estimate (Eq. 0) has to be reduced by the 
factor 3-i/("+3). 



A formula similar to Eq. Q has been derived by Brunetti 
et al. (1997V However, Eq. (A3) in Brunetti et al., which 
includes the lower energy cutoff of the cosmic ray electron 
spectrum, is not applicable in the case of dominating protons. 

3.2. Revised minimum energy formula for the total field 
From Eq. ( IA16t (Appendix) and detot / dB = we get: 



S,ni„ = { 27r(2a + 1) (a + 1) (Kq + 1) E^- 
x{iy/2cir /[(2a-l)c2(7)^C4W]} 



l/(a+3) 



(4) 



The ratio of minimum-energy magnetic and cosmic ray en- 
ergy densities is: 

ee/ecR = {a + l)/2 

Hence, the ratio of minimum-energy (|4j and equipartition (0 
field strengths is not constant, as in the classical case (Q, but 
depends on the synchrotron spectral index a: 

Bmin/Scq = [(« + l)/2]l/("+^' 

For a = 1 the revised formula gives identical results for the 
equipartition and minimum energy conditions. 

The revised formula is not valid for a < 0.5 (7 < 2) be- 
cause the second integral in Eqs. iA5\ and ( IA12t diverges 
for E2 ^ 00. Such flat injection spectra are observed in 
a few young supernova remnants and can be explained by 
Fermi acceleration operating in a strongly magnetized (low- 
f3) plasma (Schlickeiser & Fiirst 1989 1. In the diffuse medium 
of galaxies, radio lobes and clusters, a < 0.5 is observed only 
at low frequencies and indicates strong energy losses of the 
electrons due to ionization and/or Coulomb interaction (Pohl 
I1993al Sarazin fT999t where the formula cannot be used. 

Strong shocks in a non-relativistic /3 > 1 plasma may 
generate injection spectra with a — 0.5 (Schlickeiser & Fiirst 
I1989> where the total cosmic ray energy is computed accord- 
ing to Eq. (IA14t and the minimum energy formula has to be 
modified accordingly: 



Brrnn - <^ 4^(a + l)(Ko + l)/.^" 



2a 



^+ln{E2/Ep) 



X (i./2ci)" / [C2(7)^C4«]} 



1/(q+3) 



(5) 



Fig- n shows the ratio q between the revised minimum- 
energy field strength Smin (Eq- 0) and the classical value 



Bciass (Eq. ( lA20b ). using ratios of Kq = 100 and /C = 100, 
respectively. For a ~ 0.5 (7 ~ 2) the ratio is almost con- 
stant and about 2. For larger values of a, the ratio becomes a 
function of the field strength i?min- For weak fields (below a 
few /iG) and values of a between ~ 0.6 and ~ 1 the revised 
value differs insignificantly (with respect to the typical errors 
of 20% - 30%) from the classical one. 

For flat radio spectra (0.5 < a < 0.6) the classical esti- 
mate is too low because the fixed upper limit i^max = 10 GHz 
used for integrating the electron energy excludes the high- 
energy part of the cosmic ray spectrum carrying a significant 
fraction of the total energy. On the other hand, the classical 



3.4 Proton-to-electron ratio and energy losses 
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Fig. 1. Ratio between the revised minimum-energy field 
strength i?niin (Eq. 0) and the classical value i?ciass 
(Eq. ( IA20> . assuming Vynin = 10 MHz and t'max =10 GHz) 
as a function of the revised field strength Bmin and the syn- 
chrotron spectral index a. A proton-to-electron number den- 
sity ratio of Kq = 100 was adopted. 



estimate is too high for steep radio spectra (a > 0.7) if the 
field strength is > 10 /^G. Here the lower integration limit 
^'min = 10 MHz corresponds to energies < 250 MeV which 
is below the lower break energy Ep ~ 938 MeV of the proton 
spectrum (see Sect. I3.4l i. so that the total cosmic ray energy 
is overestimated. 

The ratio -Bmin/^ciass depends weakly on Kq [q oc 
-2/7^ frequency limits used for the 

classical formula. For example, q is about 10% smaller for 
a = 0.51 when increasing t^max from 10 GHz to 100 GHz. 

3.3. Relativistic electron/positron pair plasma 

Jets of radio galaxies may be composed of electrons and 
positrons (Bicknell et al.'200ri. Shock acceleration acting on 
a relativistic electron/positron pair plasma generates a power 
law in momentum which leads to a break in the energy spec- 
trum at = 511.00keV = 8.187M0-'^erg,butthisbreakis 
not observable in the radio range because low-frequency ra- 
dio spectra are flattened by ionization and/or Coulomb losses 
causing a break at Ei (Ei > E^) in the stationary energy spec- 
trum (see Sect. \3Al . The revised formulae Q or can be 
applied to an electron/positron plasma by using Kq = and 
replacing the lower energy break Ep by Ei. If, however, syn- 
chrotron or inverse Compton loss is significant beyond Egyn, 



the total cosmic ray energy as to be computed by integration 
over the observed spectrum, not according to Eq. ( IA16t . 

The ratio of the revised field strength in a relativistic 
electron/positron plasma i3min,G (Kq = 0) to that in a 
proton-dominated plasma Bmin (using Eq. and assuming 
(Ko + 1) ~Ko)is: 



Br, 



c/Bn 



= {Ep/E^)-'^°/^°'+^'> 

X (^p/Si)(2"-l)/("+3) 



(6) 



where ao is the spectral index of the injection synchrotron 
spectrum. -Bmin.c/^min varics only weakly with spectral in- 
dex a and lower energy break E^. For a wide range in ao and 
a (0.5-0.75) and in E-, (10-100 MeV), Smin.c/Smin ^ 1/3. 

3.4. Proton-to-electron ratio and energy losses 

The ratio K(E) of the proton-to-electron number densities per 
particle energy interval {np/ric) depends on the acceleration 
process, the propagation and the energy losses of protons and 
electrons. In the range of particle energies Ep < E < i?syn; 
where losses are small or affect protons and electrons in the 
same way ("thin target", see below), K = Kq is constant. 

For electromagnetic acceleration mechanisms which gen- 
erate a power law in momentum, Kq depends only on the 
transition energies from the non-relativistic to the relativistic 
regime for protons and electrons, Ep and E^, and on the par- 
ticle injection spectral index 70 (Eq. iA4i : see also Bell i 19781 



U70-l)/2 



Schlickeiser|2002| p. 472): 



Ko = f^p,o/'^c,o = 
= {Ep/E,r> 



(Ep/E, 

{Ep < E < Esyn, thin target) 



(7) 



where 70 is the spectral index of the injection cosmic ray 
spectrum. 70 is related to the shock compression ratio r via 
70 = ('^ + — 1) (see Appendix). For 70 ~ 2.2, as 

expected from acceleration in supernova remnants, we get 
Ko — 100, consistent with the local Galactic cosmic ray data 
near the sun at particle energies of a few GeV (see Appendix). 



Eqs. ( IA6I 1 and are used to compute the intrinsic ratio 
/C of total energies of protons and electrons (for negligible 
energy losses): 

IC = (Sp/Se)^'"^"'/' (8) 

Fig.|2]shows the variation of Kq and JC with the compression 
ratio r in the shock. Table^gives values of Ko and IC for var- 
ious injection processes. For strong shocks in non-relativistic 
gas (r = 4, Jo — 2) IC — Ko — 40. For steeper intrin- 
sic spectra (weaker shocks) K, decreases, while Ko increases 
(Fig.|2ji. This demonstrates why the classical formula using IC 
is of little practical use: The number density ratio Kq, in the 
particle energy range relevant for the observed synchrotron 
intensity, is not proportional to the energy ratio IC. 

In galaxy clusters secondary electrons are generated by 
hadronic interactions between cosmic ray protons and nuclei 
of the hot gas and may contribute to synchrotron emission 
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3.4 Proton-to-electron ratio and energy losses 



Table 1. Cosmic-ray proton/electron number density ratio Ko ("thin target") and the ratio /C of total energy densities for 
various injection mechanisms 

CR origin K. 



Fermi shock acceleration (strong shocks, non-relativistic gas) 40 - 100 40 - 20 

Secondary electrons 100 - 300 3> 1 

Turbulence (resonant scattering of Alfven waves) ~ 100 ~ 20 

Pair plasma 



100000 
10000 
1000 
100 - 
10 

1 - 

0.1 

0.O1 ■■ 

2 2.4 2.8 3.2 3.6 4 

compression ratio r 

Fig. 2. Variation of the intrinsic ratio of proton-to-electron 
number densities Ko (thick line) and of the intrinsic ratio 
of total energy densities K, (dotted line) with the shock com- 
pression ratio r. Energy losses are assumed to be negligible. 



(Dennison'1980"l. The spectral index 70 of the electron spec- 
trum is larger (steeper) than that of the protons 7p, depend- 
ing on the model (Pfrommer & EnBlin I2004a> . As only a 
small fraction of the proton energy can be converted into sec- 
ondary electrons, the ratio of total proton/electron energies 
is /C » 1. Assuming balance between injection and losses, 
Dennison CmU) estimated /C = 5 [9 + {B/nG)% Typical 
expected values for the number density ratio in galaxy clus- 
ters are Kq = 100 - 300 above particle energies of a few 
GeV (Pfrommer & EnBUn l200?B . 

Eq. Q is valid (i.e. Ko is constant) in the energy range 
Ep < E < Esyn if the spectral indices of the electrons 
and protons are the same. This can be achieved if energy 
losses are negligible, or if the timescale for cosmic ray es- 
cape loss is smaller than that for electron loss by nonthermal 
bremsstrahlung, which is the case for low interstellar or in- 
tergalactic gas densities ("thin target", Pohl 1992)> so that 
nonthermal bremsstrahlung is unimportant for the stationary 
electron spectrum. Energy-dependent escape steepens all cos- 
mic ray spectra in the same way, 7o = 7p = 7o + where 
— A7 ~ —0.6 is the exponent of the energy-dependent es- 
cape time tesc (see Appendix), and the revised formula for 
equipartition or minimum energy can be applied. 

If, however, the gas density is high or the escape time is 
large, the object is a "thick target" for electrons, so that non- 
thermal bremsstrahlung loss dominates the electron spectrum 
at GeV energies. The slopes of the electron and the radio syn- 
chrotron spectra are the same as for the injection spectra (70 



and ao, respectively). As the loss rate increases linearly with 
gas density n, the proton-to-electron ratio increases with gas 
density. If the proton spectrum beyond Ep is steepened by 
energy-dependent escape (tesc oc E^^'^), e.g. in galaxies, 
the ratio K(E) also depends on energy: 

K{E) (X nE-^'^ (Sp <£;< £;sy„, thick target) (9) 

The result is that the ratio K(E) varies over the whole range 
of observable energies, so that the revised formula for energy 
equipartition or minimum energy should not be applied. The 
classical formula faces a similar problem as both the syn- 
chrotron intensity and the total energy ratio are affected by 
bremsstrahlung loss, but the effect on the ratio JC of total en- 
ergies is hardly predictable. 

We will now discuss the ratio K(E) for low and high elec- 
tron energies (outside the energy range Ep < E < Esyn). In 
a proton-dominated (Ko ^ 1) cosmic ray plasma, the lower 
energy limit in Eq. Q or (|9} is set by the proton spectrum 
which flattens below Ep — 938 MeV, while the electron spec- 
trum flattens below So = 511 keV (Eq. (|a3). For E < Ep 
the proton and electron spectra are not proportional, so that 
K(E) is not constant: 

K{E) = {E/EX" iEe<E< Ep) (10) 

The revised formula should not be applied for E < Ep. In a 
10 /iG magnetic field, electrons with E < Ep are observed at 
frequencies below about 140 MHz (Eq. (|2ji). 

The classical formula is also affected if radio observations 
at low frequencies are used. Even in a weak 1 /iG magnetic 
field, the minimum allowed frequency corresponding to Ep is 
still higher than the standard lower frequency limit of 10 MHz 
used in the classical formula (Eq. Q). The standard integra- 
tion range of > 10 MHz would trace an appropriate part of the 
energy spectrum only if the field strength is below 0.7 fiG. 
However, synchrotron emission from such weak fields is too 
faint to be detected with present-day radio telescopes. Hence, 
the lower frequency limit of 10 MHz used in the classical for- 
mula (Eq. Q) is generally too low and leads to an overesti- 
mate of the field strength. This bias increases with increasing 
field strength (see Fig.nj. The ratio IC of total energies cannot 
be applied to low-frequency radio data. 

At low energies, energy losses modify the proton and 
electron spectra (Pohl ll99'3al and the effective K(E) is differ- 
ent from that according to Eq. dlOt . The break in the electron 
spectrum by A7C = — 1 (flattening) due to ionization and/or 
Coulomb loss depends on gas density (see Eq. (1) in Brunetti 
et al.' 19971 and occurs at a few 100 MeV in the ISM of galax- 
ies and at a few tens MeV in radio lobes and clusters. These 
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energies correspond to radio frequencies which are generally 
below the observable range, so that this effect is irrelevant 
for equipartition estimates in a proton-dominated relativistic 
plasma. 

E'syn in Eq. (0 is the upper energy break where syn- 
chrotron or inverse Compton loss of the electrons becomes 
dominant. Inverse Compton loss has the same energy and fre- 
quency dependence as synchrotron loss; both are of similar 
importance if the energy density of the radiation field and the 
magnetic energy density are comparable. Inverse Compton 
loss from the CMB background dominates if the field strength 
is below 3.25 x (1 + z)^ ^G. For Ec > i?syn the spectral in- 
dex steepens by A70 = 1, observable as a smooth steepening 
of the synchrotron spectrum by Aa = 0.5. In the revised for- 
mula for energy equipartition or minimum energy one should 
not use data at radio frequencies corresponding to the elec- 
tron energy range E > Egyn, where their energy spectrum is 
not proportional to the proton spectrum and the ratio K is a 
function increasing with E (Pohl ll993al Lieu et al. lT999l l. A 
simplified estimate is: 

KiE) =Ko{E/Esyn) {E>Esyn) (11) 

Using instead Kq in the revised formula leads to an underes- 
timate of the field strength. Again, the classical formula has a 
similar problem because the ratio K, of total energies also in- 
creases in the case of strong synchrotron or inverse Compton 
loss. 

If the field strength varies along the line of sight or within 
the volume observed by the telescope beam, synchrotron loss 
may lead to an anticorrelation between field strength and cos- 
mic ray electron density, so that the equipartition field is un- 
derestimated further (Beck et al. .2003j . 

Synchrotron loss is significant in sources with strong 
magnetic fields like young supernova remnants, starburst 
galaxies and radio lobes (Sects. 14.21 and 14. 3> . and also in 
galaxies away from the acceleration sites of cosmic rays, for 
example in interarm regions, in the outer disk and in the halo 
(Sect.im. 

The various bias effects are summarized in Table |2] 



3.5. Equipartition and minimum-energy estimates of the 
regular field 

Knowing the equipartition or minimum-energy estimate of 
the total field strength, the equipartition or minimum-energy 
estimate of the regular field strength -Biog,± in the sky plane 
can be computed from the observed degree of polarization p 
of the synchrotron emission: 

7 1 

P - Po (1 + ^q') I (1 + + -let) (12) 

where po is the intrinsic degree of polarization (po = (3 — 
3a)/ (5 — 3a)) and q is the ratio of the isotropic turbulent field 
Stuib to the regular field component -Biog,_L in the sky plane 
(Beck et al. l2003t . For the case of a dominating turbulent field 
(g> l),p~ 2.1po(?"^. 



4. Discussion and examples 
4.1. Weak fields: spiral galaxies 

The Milky Way and most spiral galaxies have steep radio 
spectra in the frequency range of a few GHz which indicates 
that the energy spectra of their cosmic ray electrons are steep- 
ened by escape loss (see Appendix). Galaxies are "thin tar- 
gets" for cosmic ray electrons, except for massive spiral arms 
and starburst regions (see Sect. 14. 2t . The revised formula can 
be applied, using the part of the radio spectrum below the 
break frequency i/sy„ beyond which synchrotron or inverse 
Compton losses become important. The upper energy limit 
i^syn corresponding to i/sy„ can be estimated as follows. 

The synchrotron lifetime igyn for electrons of energy E is 
(Lang[T953 p. 32): 

i,yn = 8.35-109yr/[(Bi^G)2 (£;,y„/GeV)] 

= 1.06-109yr/[(Bi/^G)i-5 (;/syn/GHz)0-5] (^3) 

Note that the constant is different from that used in most pa- 
pers (e.g. Carilli et al. [79911 Mack et al. [T998l Feretti et al. 
1998 1. The escape time tcsc in the Milky Way is ~ 10'' yr at 
non-relativistic and mildly relativistic energies (Engelmann 
et al. 1990; Schlickeiser 2002, p. 439), so that the break in 
the radio frequency spectrum is expected at: 

i/,y„ « 10GHz/(Si/10AiG)3 (14) 

For the typical strength of ISM magnetic fields of 10 pQ 
(Niklas 1995) we get i/gyn ~ 10 GHz. At higher frequencies 
the synchrotron intensity ly and the spectral index a can be 
used neither for the revised formula nor for the classical for- 
mula. If the resulting field strength is larger than 10 pQ, like 
in grand-design spiral galaxies, the useful frequency range 
shifts to even lower frequencies. Furthermore, the thermal 
contribution to the total radio intensity increases with increas- 
ing frequency and has to be subtracted. 

In galaxies with a high star-formation rate, galactic winds 
may form. The escape time may then be shorter than in the 
Milky Way and j/gyn is larger than according to Eq. ( I14t . Sig- 
nature of galactic winds is that at low frequencies the local 
radio spectrum in edge-on galaxies remains flat up to large 
heights above the disk (Lerche & Schlickeiser 1982). In the 
Milky Way, the radio spectral index distribution is consistent 
with a galactic wind (Reich & Reich fl988l l. 

Pohl ( 1993ai argued that in the spiral galaxy M 51 the 
spectral index of the total radio emission (integrated over the 
whole galaxy) is a ~ 0.5 (70 — 2.0) at low frequencies, but a 
frequency break at about 1 .4 GHz due to energy losses causes 
a steepening of almost the whole observable radio spectrum. 
As a result, the classical minimum-energy field strength in 
the spiral galaxy M 5 1 is too high. However, the compilation 
of all data by Niklas J1995t does not indicate a significant 
steepening of the radio spectrum of M 5 1 until 25 GHz, but 
a flattening due to thermal emission. In his sample of spiral 
galaxies, Niklas J1995> found only very few cases of spec- 
tral steepening. Hence, escape seems to be the dominant loss 
process in the spectra of integrated radio emission of galaxies 
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4.2 Strong fields: young supernova remnants and starburst galaxies 



Table 2. Bias of equipartition field estimates. K(E) is the cosmic ray proton/electron ratio, where Kq denotes a value which 
does not vary with energy 



Method 


Cosmic-ray composition 


Bias effect 


Field strength 


Reference 


Classical 


+ e (Ko ~ 100) 


Integration over frequency 


Underestimate (a < 0.6) 
Overestimate (a > 0.7) 


This paper 
(Fig.El 


Classical 


e" + e+ (Ko = 0) 


Fixed frequency range 
(e.g. 10 MHz-10 GHz) 


Underestimate (weak fields) 
Overestimate (strong fields) 


This paper 


Classical+revised 


p+ + e- (K(E)> 100) 


Synchrotron/IC/ 
bremsstrahlung losses 


Underestimate 


This paper 


Classical+revised 


any 


Field fluctuations without/ 
with synchrotron loss 


Overestimate (weak fields)/ 
Underestimate (strong fields) 


Beck et al.<2003l 



between about 100 MHz and 10 GHz. (This does not hold for 
spectra of the local emission, see below.) 

Fitt & Alexander ( 1993 1 used the 1.49 GHz radio fluxes 
of 146 spiral and irregular galaxies to derive minimum- 
energy field strengths with the classical formula, simply as- 
suming a constant spectral index of a = 0.75, a negligible 
thermal fraction and the standard frequency limits. Their dis- 
tribution peaks around 1 1 jiG with a standard deviation a of 
4 fjG (for a ratio of total energies JC of 100). 

Niklas J 19951) observed a sample of spiral galaxies in ra- 
dio continuum at 10.55 GHz and compiled all available data 
at other frequencies. Based on the spectra of the integrated 
radio emission, Niklas et al. (|T997i could separate the syn- 
chrotron from the thermal emission for 74 galaxies. The mean 
thermal fraction at 1 GHz is 8%, and the mean synchrotron 
spectral index is a = 0.83 with a standard deviation of 0.13. 
Hence, the average spectrum of the extragalactic cosmic ray 
electrons (and probably also that of the protons) at particle 
energies of a few GeV has the same spectral index (7 ~ 2.7) 
as that in the Milky Way. Niklas (Il995t derived galaxy- 
averaged field strengths according to the classical formula 
(IA20> and to a revised formula similar to Eq. 0, approximat- 
ing integration (IA12> by one integral from Ei = 300 MeV 
to infinity. His results (for Kq = 100) are similar for the two 
cases, as expected from Fig. [2 The mean is 9 fiG (with a 
standard deviation a — "i fiG) both for the revised formula 
and for the classical case. The distribution of field strengths 
is more symmetrical for the revised case. 

However, application of the equipartition formula to re- 
gions within spiral galaxies or in their radio halos needs spe- 
cial care. For example, the "ring" of M 31 emitting in radio 
continuum is formed by cosmic ray electrons, while the mag- 
netic field extends to regions inside and outside the ring, as 
indicated by Faraday rotation data (Han et al. I1998> . Leav- 
ing the star-formation regions, cosmic ray electrons rapidly 
lose their energy due to synchrotron loss. Hence, K(E)> 100, 
so that the equipartition formula underestimates the field 
strength inside and outside the ring. The same argument holds 
for the outer disk and the halo of galaxies, far away from the 
acceleration sites of cosmic rays. 

In galaxies with strong magnetic fields like M 51, syn- 
chrotron loss may dominate already at a few 100 pc from the 
acceleration sites, e.g. in interarm regions, as indicated by 



spectral steepening between the spiral arms (Fletcher et al. 
I2005> . Hence, the interarm field strengths given by Fletcher 
et al. f2004i) are underestimates. 

4.2. Strong fields: young supernova remnants and 
starburst galaxies 

The classical equipartition formula overestimates the mag- 
netic field strength in objects with proton-dominated cos- 
mic rays and with strong fields (Fig. [0, like massive spiral 
arms, young supernova remnants or starburst galaxies. Here 
the lower frequency limit is too low and hence overestimates 
the total cosmic ray energy (see Sect. 13. 4> . In a 100 /iG field, 
for example, the energy of an electron radiating at 10 MHz is 
only 80 MeV, where the corresponding proton number in any 
medium is strongly reduced (E ^ Ep) and the proton-to- 
electron number density ratio K(E) may drop to values even 
smaller than 1 . 

Hence, many values for the field strength quoted in the lit- 
erature are strong overestimates. For example, Allen & Kro- 
nberg ( ,1998j integrated the radio spectrum from 10 MHz to 
100 GHz of young supernova remnants in the galaxy M82 
and derived field strengths of a few mG for /C = 0, scal- 
ing with (/C + 1)^/^. The revised formula gives values which 
are between 4 and 50 x smaller, depending on spectral index 
a. Perez-Torres et al. J2002> derived a minimum-energy field 
strength in SN 1986J (a = 0.69) of about 13 mGx (/C + 1)2A 
which, according to our revised formula, is too high by about 
one order of magnitude. On the other hand, the equipartition 
strength of the flat-spectrum SN 1993J (a = 0.51) of 38 mG 
(Chandra et al. 20Q4} is too low; the revised formula gives 
~ 100 mGx (JC + Xfl'^ which fits better to the field strength 
derived from the synchrotron break energy (by inserting the 
SN age as tsy„ in Eq. (fT3l). 

Note that the cosmic ray electron distribution in young 
supernova remnants is not in a stationary state, so that energy 
losses modify the spectrum in a different way than discussed 
in Sect. 13. 41 The observed spectrum is that of injection, possi- 
bly modified by synchrotron loss, so that the revised formula 
can be applied. 

Hunt et al. (I2004> discussed the radio emission of the blue 
compact dwarf galaxy SBS 0335-052. Applying the classi- 
cal formula (assuming = 40 and a frequency range of 



4.3 Radio lobes and galaxy clusters 
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10 MHz-100 GHz) these authors obtained an equipartition 
field of ~ 0.8 mG, while the value from our revised for- 
mula is about 30 /iG which fits much better to other starburst 
galaxies hke M82 (Klein et al.Q,988 1 or the Antennae (Chyzy 
& Beck 1200 4^ or blue compact dwarf galaxies (Deeg et al. 
119931 . The same discrepancy may arise in galactic nuclei 
with starbursts. Beck et al. (I2005> derived field strengths of 
~ 60 in the central starburst regions of the barred galax- 
ies NGC 1097 and NGC 1365, while the classical formula 
would give much larger values. 

However, starburst galaxies and regions of high star for- 
mation rate in the central regions and massive spiral arms of 
galaxies have high gas densities and mostly flat radio spectra 
(a = 0.4-0.7), so that they probably are "thick targets" for 
the cosmic ray electrons. If so, the equipartition estimate is 
too low, and the correct value can be computed only by con- 
structing a model of gas density and cosmic ray propagation. 

4.3. Radio lobes and galaxy clusters 

The classical equipartition field strengths are 20-100 fiG in 
the radio lobes of strong FRJl radio galaxies (Brunetti et al. 
EW), 100-600 fiG in hot spots (Carifli et al. 1991 Hams et 
al. i2000. Wilson et al..2000. Wilson et al. .2001t Brune tti et 
al. l2001al Hardcastle l200Tl Hardcastle et al. l200ni2002t . and 
0.5-1 /iG in the diffuse gas of galaxy clusters and relics (Fer- 
etti & Giovannini 1996j . In all these cases it was assumed that 
relativistic protons are absent (K, ~ 0, "absolute minimum 
energy") or contribute similarly to the total energy density as 
the electrons (K. — V). 

According to Fig. |2] /C < 1 could also be the result of 
a weak shock. However, particle acceleration at low Mach 
numbers is inefficient (Bogdan & Volk <1983i Volk et al. 
I1988I I. so that its contribution to the total cosmic ray popu- 
lation is negligible. 

If /C = 0, the radio spectrum reflects the energy spec- 
trum of the total cosmic rays, and the equipartition formula 
gives reliable results, supposing that the integration limits are 
set properly (Sect. l33l . Integration over the range of Lorentz 
factors of the electrons/positrons according to the breaks as 
observed in the radio spectrum of sources leads, in the case 
of weak magnetic fields, to slightly larger field strengths than 
applying the classical formula with fixed frequency limits 
(Brunetti et al. ll997t . However, the lower break energy Ei in 
the electron/positron energy spectrum needed for the equipar- 
tition estimate may correspond to a radio frequency which is 
too low to be observed. For example, electrons with the min- 
imum Lorentz factor of 50 assumed for the lobes of 3C219 
(Brunetti et al. 1999 1 in a 10 ^G field radiate at 100 kHz, 
a frequency much below the observable range. Here the fre- 
quency limit of z^min = 10 MHz assumed in the classical 
method is too high, and hence the field strength is underesti- 
mated, by 1.5 X and 2x for a — 0.8 and a = 1, respectively, 
or by an even larger factor if the minimum Lorentz factor is 
smaller than 50. 

If the fields are strong, the break energy Ei of the electron 
spectrum is observable in the radio spectrum. Hardcastle et al. 
(12002^ determined typical minimum Lorentz factors of 1000 



from the breaks observed at 0.5-1 GHz in hot spots of radio 
lobes and estimated field strengths of 100-200 ^G (assum- 
ing equipartition with Kq = 0). In such cases, the classical 
formula gives overestimates because i^niin = 10 MHz is too 
small. 

The occurence of electron/positron plasmas {K, — 0) 
in astronomical objects is under debate. A relativistic elec- 
tron/positron plasma was discussed for the jets (Bicknell et 
al. 2001 1. However, observational evidence tends to favour 
an electron/proton plasma (Leahy et al. 1200 U . Even if jets 
do eject electrons and positrons into a cluster with the same 
energy spectrum, these particles will not survive long enough 
to generate significant radio emission from the IGM medium, 
so that (re)acceleration by intergalactic shocks (Blandford 
& Ostriker 1978, Sarazin 1999, Gabici & Blasi 200 3), by 
interaction with MHD waves (Schhckeiser & Miller [19981 
Brunetti et al. '2001b' Fujita et al. '2003 1, or by reconnection 
(Hanasz & Lesch 2003 1, or production of secondary electrons 
(Pfrommer & EnBlin l2004al IS necessary. 

If, however, the contribution of relativistic protons to 
the total cosmic ray energy is dominant (/C ^ 1), as pre- 
dicted by all electromagnetic acceleration models (Table [0, 
the equipartition field strength increases with respect to that 
for /C = 0. In the hot spots of radio lobes (Wilson et al. 
l200n Hardcastlel200n Hardcastle et al. 2001 ,,.200 2 1 stronger 
equipartition fields would fit much better to the high field 
strengths derived from Faraday rotation measures (Feretti et 
al. [T995l Fusco-Femiano et al. 1999.1 . 

Another reason for systematic field underestimates are 
energy losses. The electron spectrum is a result of the bal- 
ance between acceleration and energy losses. Ionization loss 
dominates at low energies; at high energies synchrotron loss 
dominates in clusters with magnetic fields stronger than a few 
fxG while inverse Compton loss due to CMB photons domi- 
nates in clusters with weaker fields (Sarazin 1999V E > Egyn 
is the energy range where synchrotron and/or inverse Comp- 
ton losses are strong and the electron spectrum is steeper than 
the proton spectrum (70 = 7p + 1). The proton-to-electron 
number density ratio K increases with energy according to 
Eq. il 1\ . In this case, the revised formula underestimates the 
field strength. For example, using a value for K of 100 instead 
of 1000 would underestimate the equipartition field strength 
by a factor of about 2. The classical formula has a similar 
problem. 

Synchrotron and inverse Compton losses steepen the ra- 
dio spectra. The radio spectra of most radio lobes and hot 
spots in the range 0.4-15 GHz are well approximated by 
power laws with spectral indices between of a ~ 0.7-0.8 
in the frequency range 0.4-15 GHz, with a spectral break to 
a > 1 between the radio and X-ray regimes. Meisenheimer 
et al. ( 1989 1 fitted the spectra of five bright hot spots with a 
low-frequency spectral index of a ~ 0.5 with smooth breaks 
which start to steepen the spectra already beyond a few GHz. 
The spectral index of the diffuse radio emission from the 
Coma cluster is a ~ 1.0 below 100 MHz and strongly steep- 
ens at higher frequencies (Fusco-Femiano et al. 119991 Thier- 
bach et al. .2003j . which is a clear sign of energy loss. The 
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4.4 Comparison with field estimates .. 



(revised) equipartition field strength in the Coma cluster is 
0.6 /zGx (K + 1) V("+3) (Thierbach et al.'^OOS^, which gives 
~ 4 for a = 0.8 and K = 1000. This value is not far 
from that derived from Faraday rotation data (Feretti et al. 
[T995t . 

Pfrommer & EnBlin J2004bt modeled the radio spectrum 
of the Coma and Perseus clusters assuming equilibrium be- 
tween injection of secondary electrons (hadronic interaction) 
and energy losses. Their self-consistent minimum-energy for- 
mula, which does not require an assumption about K. or K, 
gives slightly larger field strengths for the Coma and Perseus 
clusters compared with the classical formula. However, this 
model cannot account for the radial steepening of the radio 
spectra observed in both clusters. 

Our revision of total field strengths may have significant 
effects on estimates of the age of electron populations based 
on the synchrotron break frequency in the radio spectrum 
(Eq. il3\ ). Application of the classical formula for objects 
with strong fields may significantly overestimate the field 
(Fig-E ™d hence underestimate the electron lifetime. More 
serious is the assumption of a proton/electron total energy ra- 
tio = 1 in most papers (Carilli et al. '1991' Mack et al. 
119981 Feretti et al. 1998 Parma et al. 1999 1 which is not sup- 
ported by any acceleration model (see above). Kq ~ 100 
in the revised formula would increase field strength by 3.5 x 
and hence decrease the synchrotron age by 6.5 x. The dis- 
crepancy between dynamical and synchrotron ages found by 
Parma et al. (I1999t is increased. This problem needs to be 
re-investigated. 

4.4. Comparison with field estimates from tlie ratio of 
synclirotron to inverse Compton X-ray intensities 

The same cosmic ray electrons emitting radio synchrotron 
emission also produce X-ray emission via the inverse Comp- 
ton effect. The ratio of radio to X-ray intensities can be used 
as a measure of field strength (Harris & Grindlav ll979l Gov- 
oni & Feretti l2004> . Comparison with the equipartition esti- 
mates reveals three cases: 

(1) A high equipartition strength is needed to achieve a 
consistent picture. 

The lobes of two low-power radio galaxies reveal an ap- 
parent deficit in X-ray emission compared with radio emis- 
sion when assuming no protons (JC = 0) (Croston et al. 
|2^3j . Relativistic protons with 300-500 times more en- 
ergy than the electrons would increase the equipartition field 
strength from ~ 3 fiG to ~ 15 /iG and reduce the number 
density of electrons, which would explain the weak X-ray 
emission and also ensure pressure balance between the radio 
lobes and the external medium. 

(2) A low equipartition strength is needed to achieve a 
consistent picture. 

The equipartition values (assuming JC ~ or 1) are 
similar to those derived from X-ray emission by inverse 
Compton scattering of the background photons in most ra- 
dio lobes (Brunetti et al. ■1997.) and in the Coma cluster 



(Fusco-Femiano et al. 119991 see also Table 3 in Govoni & 
Feretti 2004). In most hot spots the equipartition values are 
also similar to those derived from X-ray emission by inverse 
Compton scattering of the synchrotron photons ("synchrotron 
self-Compton emission"), e.g. in Cyg A (Harris et al. 119941 
Wilson et a l. 2000 1, 3C123 (Hardcastle et_al. 200D, 3C196 
(Hardcastle l2noll . 3C295 ( Harris et al. l2000t . in 3C263 and 
in 3C330 (Hardcastle et al. '2002't. If, however, /C > 1, as 
prediced by acceleration models (Tabled, the field estimates 
increase and fit better to the Faraday rotation data (see above), 
but contradict the X-ray data in these objects which then are 
"too bright" in X-rays. 

(3) The equipartition strength is always too high. 

Some objects are "too bright" in X-rays compared with 
the energy density of the cosmic ray electrons derived from 
the equipartition assumption with /C = 0. In other words, the 
equipartition values even for /C = are already several times 
larger than those allowed from the ratio of radio/X-ray inten- 
sities. This is the case in the lobes of PKS 1343-601 (Tashiro 
et al. a998j and 3C219 (Brunetti et al. .l999 l. and in the hot 
spots of Pictor A (Wilson et al. 200^, 3C351 (Brunetti et al. 
2001ai and 3C351 (Hardcastle et al. l2002> . Dominating pro- 
tons (/C ^ 1) would even increase the discrepancy. 

The correction introduced for strong fields by our revised 
equipartition formula (Fig.[0 cannot solve the problem. For 
example, in the bright core of the western hot spot of Pictor 
A, Wilson et al. (2001 1 found the largest discrepancy (more 
than 10) between the field strength derived from the classi- 
cal equipartition formula (~ 470 /iG, assuming IC ~ 1) and 
that from the radio/X-ray intensity ratio (~ 33 fiG). Our re- 
vised formula (for Kq = 1) reduces the equipartition value 
only to ~ 350 /iG, not enough to remove the discrepancy. If 
cosmic ray protons dominate, the equipartition field strength 
increases to ~ 1 mG (for Kq = 100 and a = 0.8), so that 
the discrepancy between the field estimates increases further. 

Carilli & Taylor (I2002> discussed possible solutions. 
Firstly, most of the X-ray emission may not be of inverse 
Compton, but of thermal origin, which should be tested with 
further observations. Secondly, an anisotropic pitch angle dis- 
tribution of the cosmic ray electrons could weaken the syn- 
chrotron relative to the inverse Compton emission. Thirdly, 
in magnetic fields of about 1 /iG the observed synchrotron 
spectrum traces electrons of larger energies than those emit- 
ting the inverse Compton spectra. If the electron spectrum 
steepens with energy (e.g. due to energy losses), the radio 
emission is reduced at high energies. 

Another possible (though improbable) explanation of the 
discrepancy is that equipartition between magnetic fields and 
cosmic rays is violated. An increase in electron density by 
typically n'^/uc^cq — 5 is required to match the excessive X- 
ray intensities in the objects listed above. The corresponding 
decrease in field strength for a fixed radio synchrotron inten- 
sity and a ~ 0.8 is B'/B^q = «/nc,oq)"^/'"+^' ^ 0.4, 
and the energy density ratio q between cosmic rays and mag- 
netic fields increases by q = (rio/riccq)^"^'^''^'"^^-' — 30. 
Such an imbalance between particle and field energies is un- 
stable and would cause a rapid outflow of cosmic rays. 
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Finally, the magnetic field may be concentrated in fil- 
aments. Intracluster magnetic fields can be amplified by 
shocks in merging clusters (Roettiger et al. 1999 1. The ra- 
tio of synchrotron to X-ray intensities from the same cosmic 
ray electron spectrum is Isyn/Ix oc < n^B^'^"' > / < rto >• 
If small-scale fluctuations in and in B are uncorrelated, 
Isyn/Ix oc < > = < B >^ //b, where /b is the effec- 
tive volume filling factor of the magnetic field, so that the 
field strength estimate from the radio/X-ray intensity ratio 
varies with /~°'^. To explain a lOx discrepancy between the 
field estimates, a filling factor of 10^^ is required. The en- 
hanced magnetic field in the filaments leads to synchrotron 
loss of the electrons and hence to an anticorrelation between 
rtc and B. The X-ray emission may be biased towards regions 
of weak fields, while synchrotron intensity is biased towards 
regions of strong fields. This explanation appears plausible 
and deserves further investigation. 

5. Summary 

Assuming that cosmic rays are generated by electromagnetic 
shock acceleration producing the same total numbers of rela- 
tivistic protons and electrons with power laws in momentum, 
we showed that the ratio of total energies of protons and elec- 
trons K, is about 40 for strong shocks and decreases with de- 
creasing shock strength. The ratio K of number densities of 
protons and electrons per energy interval for particle energies 
E >1 GeV is also about 40 for strong shocks and increased 
with decreasing shock strength. Both ratios further depend on 
the various energy loss processes for protons and electrons. 

The classical equipartition or minimum-energy estimate 
of total magnetic fields strength from radio synchrotron in- 
tensities is based on the ratio /C of the total energies of pro- 
tons and electrons, a number which is hardly known because 
the proton and electron spectra have not been measured over 
the whole energy range. We propose a revised formula which 
instead uses the number density ratio ratio K in the energy 
interval relevant for synchrotron emission. This ratio can be 
adopted from the value observed in the local Galactic cosmic 
rays, or it can be estimated for the energy range where energy 
losses are negligible or particle escape loss dominates ("thin 
target"), so that K is constant (Kg). 

Furthermore, the classical equipartition / minimum- 
energy estimate is incorrect because the cosmic ray energy 
density is determined by integration over a fixed interval in 
the radio spectrum, which introduces an implicit dependence 
on the field strength. Applying our revised formula, the field 
strengths for a proton-dominated relativistic plasma are larger 
by up to 2 X for flat radio spectra (synchrotron spectral index 
a < 0.6), but smaller by up to lOx for steep radio spectra 
(a > 0.7) and for total field strengths S > 10 /iG (Fig.[T]and 
Table|2j- The revised formula can be applied if energy losses 
are negligible or if escape is the dominant loss process for 
cosmic ray electrons at GeV energies. The classical field es- 
timates for young supernova remnants are probably too large 
by a factor of several. The average field strengths for spi- 
ral galaxies remain almost unchanged. For objects containing 
dense gas ("thick targets") where energy loss of cosmic ray 



electrons by nonthermal bremsstrahlung is significant, e.g. in 
starburst galaxies, massive spiral arms or in starburst regions 
in galaxies, neither the revised nor the classical estimate can 
be applied. 

Equipartition values for radio lobes and galaxy clusters 
are usually computed assuming a constant cosmic ray pro- 
ton/electron energy ratio of /C = (i.e. pure electron/positron 
plasma) ot K, — 1 (i.e. the same contribution of protons and 
electrons to the total cosmic ray energy). However, all current 
models of cosmic ray origin predict that protons are dominant 
(/C > 1, Ko ~ 100-300), so that the field estimate is too low 
by a factor (Kq + I)^^'-"^'^''. Furthermore, the radio spectra of 
radio lobes and clusters indicate that synchrotron or inverse 
Compton loss of the cosmic ray electrons are significant, so 
that K, and K and hence the field strength increase further 
The revised, stronger fields in clusters are consistent with the 
results of Faraday rotation observations. 

In case of strong fields in radio lobes, the discrepancy 
with the much lower field estimates from the radio/X-ray in- 
tensity ratio in several hot spots cannot be solved and requires 
alternative explanations, e.g. a concentration of the field in 
small filaments with a low volume filling factor, or a thermal 
origin of the X-ray emission. 

Our code B FIELD to compute the strengths of the total 
and regular fields from total and polarized synchrotron inten- 
sities, allowing various assumptions as discussed in this pa- 
per, is available from www.mpifr-bonn.mpg.de/stajf/mkrause. 
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Appendix A: The revised formula for energy 
equipartition or minimum energy 

In all electromagnetic acceleration models, like the first-order 
and second-order Fermi shock acceleration (Drury 1983J 
Schlickeiser l2002l Gabici & Blasi'2003'l and the acceleration 
by MHD waves (Schlickeiser & Miller 1998 Brunetti et al. 
2001b i Fujita et al.'2003 1, the energy spectrum of cosmic rays 
is a power law in momentum which translates into a power 
law in particle energy E with a break at E^, = mc^, where m 
is the rest mass of the accelerated particle (Bell ll978t : 

n{E)dE = no {E^,/ Eq)'-^-^"^/'^ {E/ Eo)-'^^"^^^/'^ dE 

{Ei<E< Sb) , (Al) 

n{E)dE = 710 {E/Eo)-'^°dE 

(E > E^) (A2) 
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n{E) is the number of cosmic ray particles per unit volume 
and per unit energy interval. i?o {Eq > Eh) is an energy nor- 
malization to obtain correct units; uq — n{Eo); Eo=l erg if 
all calculations are performed in cgs units. 70 is the injection 
spectral index of the energy spectrum. Ei is the threshold en- 
ergy for acceleration (see below). 

The solution of the transport equations for Fermi acceler- 
ation in a strong, non-magnetic shock in a non-relativistic gas 
(compression ratio r = 4) yields an injection spectral index 
of 7o = (r + 2)/(r - 1) = 2.0 and 70 = 2.5 in a rela- 
tivistic gas (compression ratio r = 3) (Drury[l983j Schlick- 
eiser 2002 1. The limit for ultrarelativistic shocks is 70 = 2.23 
(Kirk et al. (2000j). In a low-/3 plasma (i.e a dominant mag- 
netic field, where (3 is the ratio of thermal to magnetic energy 
densities) the compression is larger and 70 can be flatter than 



For El ^ Eh and 2 < 70 < 3, the high-energy particles 
carry the total energy and Eq. iA5i gives: 



2 (Schlickeiser & Furst ll989l l. If P is close to unity, the in- 
jection spectral index is 70 ~ 2, almost independent of the 
compression ratio. 

In models of first-order Fermi acceleration in supernova 
shocks the effective injection spectral index of protons (and 
probably also electrons) is 70 ~ 2.1-2.3 (Bogdan & Volk 
119831 Volk et al. I1988> which agrees well with the source 
spectrum fitted to the observations of local Galactic cosmic 
ray nuclei near the sun (Engelmann et al. 1990 1. Intergalactic 
shocks due to merging clusters produce electron spectra with 
70 ~ 2.3-2.4 (Gabici & Blasi'^031. 

The total number of cosmic ray particles is (for 7 > 1 and 
negligible energy losses): 

N ^ n{E)dE 

J El 



70-1 
2noEo 



[2{EhEi/El)^^-^o)/2 _ (i5^/i5o)(i-7o)- 



{El/EhEi)'^^'>-^y^ 



(A3) 



70-1 

where the approximation is valid for Ei <C ii'b- 

We adopt the commonly used assumption that the same 
total numbers of non-relativistic protons and electrons are ac- 
celerated above a certain energy threshold Ei, e.g. the non- 
relativistic kinetic energy of i?i ~ 10 KeV. From this require- 
ment it follows that the ratio of number densities at energy Eq 
(Eq > Ep) has to be: 

np,o/ne.o = (Ep/E,)^''"-^^/^ (A4) 
where Ep and Eg are the spectral break energies for protons 
and electrons, respectively. 

(Ep = 938.28 MeV = 1.5033 • lO^^ erg and E^ = 511 keV 
= 8.187- 10-^ erg) 

The total energy density of cosmic ray particles (for neg- 
ligible energy losses) is: 



e = / n{E)EdE 
Iei 

no {Eh/E,)^^-^o)/2 (i5/i5„)(i-7o)/2 

IE, 



nOO 

/ noEo{E/Eo)^-^"dE 

JEt, 



noEliEo/Eh)^'>~^{: 



1 



(A6) 



3-70 70-2 
For 7o = 2 (very strong shocks): 

e - no -Bo' ( 2 + HE2/Eh) ) (Al) 

where E2 is the high-energy limit of the spectrum according 
to diffusion coefficient and/or the finite age of the accelera- 
tion region (Lagage & Cesarsky 1983j. The TeV 7-ray image 
of a supernova remnant confirmed the existence of particles 
of very high energies (Aharonian et al. l2004> . For 70 = 3: 

e = no El {Eo/Eh) ( 1 + HEh/E^) ) (A8) 

and for 70 > 3 (weak shocks) the low-energy particles carry 
the total energy: 

e = n„El{Eo/Eh)^''°-'^/'{Eo/Ei)('">-'^/^ 

X (t^^) (A9) 

3-70 

As a consequence of the different rest masses of pro- 
tons and electrons, their spectral break energy Eh differs, and 
hence the energy spectrum of electrons steepens at lower en- 
ergies than that of protons (see Fig. 1 in Pohl I993a). The 
proton number per energy interval dominates for E > E^ = 
511 keV (see Eqs. Q and ([lO)). 

Eq. iA6\ is valid only for the injection spectrum. The ob- 
served spectral index beyond 10 GeV is 7 ~ 2.7 for many 
cosmic ray species near the sun (Fulks I1975> and in spiral 
galaxies (Niklas et al. I1997> . which is steeper than the in- 
jection spectrum and indicates an energy-dependent escape 
time from the Milky Way of ^csc oc E^^-^ (Engelmann et 
al. I1990> . Furthermore, the sub-relativistic cosmic ray pro- 
tons suffer from ionization and/or Coulomb losses, so that 
the spectrum becomes almost flat at energies below Ep with 
an energy spectral index of (2 — 7p o)/2 ~ 0.1-0.2 (e.g. Pohl 
'1993al. The observed spectrum of local Galactic protons, cor- 
rected for modulation by the solar wind, gradually steepens 
betwe en 100 MeV and 10 GeV (Fulks [D751 Ip & Axford 
I1985> and can be described as: 

n{E)dE = const dE = np,o (Ep/Eo)-'^^ dE 

(i;<100MeV) (AlO) 

n{E)dE = Hpfl {E / Eq)-''^ dE 

(-B^ 10 GeV) (All) 

with 7p ~ 2.7. A similar spectral behaviour of the protons 
is expected in all objects containing ionized or neutral gas, 
so that Eq. (|^n) is of general applicability. Replacing the 
observed gradual steepening around Ep by a break at Ep, the 
energy density of the protons is approximated by: 



en ■= / npfl{Ep/Eo)-^- EdE 

I El 



/•oc 

(A5) + / npflEo{E/Eo)^-'^-^dE 

Je„ 



(A12) 
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For El < Ep and 7p > 2: 



2 7p-2 



(A13) 



For very strong shocks which generate injection spectra 
with 7p — 2.0 (Sect. ll!2t . Eq. ( IA13t has to be replaced by: 



Ep = "P^o (2 + HE2/Ep) ) 



(A14) 



The scaling parameter ripfi of the proton spectrum can be 
determined in "thin targets" from that of the electron spec- 
trum riefi and the proton-to-electron number density ratio 
Ko, and ne,o can be replaced by the synchrotron intensity 
(surface brightness) (in erg cm~^ Hz~^ sterad"^), 
observed at the frequency i/, and the magnetic field strength 
B± in the sky plane (in G), making use of the synchrotron 
formula (e.g. Pacholcz vk l 1 970t : 



I. - C2(7c) ne,o E^^ I (i./2ci)(i->)/2 ^ 



(7c + l)/2 



(A15) 



where 

Cl 

C2(7o) 
C3 



3e/(47rmo^c'') = 6.26428 • 10'" erg-^ g-i 



(7c + 7/3)/(7e 



l)r[(37c-l)/12] 
X r[(37c + 7)/12] 
= \/3eV(47rmcc2) 
= 1.86558 • 10"^^ erg G~^ sterad"' 



C2(7c) is identical to 05(70) in Pacholczyk J 19701 p. 95 and 
tabulated on p. 232). nc{E) is the number density of cosmic 
ray electrons per unit energy interval (in cm^'^ erg^^) in the 



relevant energy range. (Note that Pacholczyk J1970> used a 
different definition of ndE) which leads to strange units.) I 
is the pathlength through the emitting medium. The spectral 
index a of the synchrotron emission relates to the spectral 
index of the energy spectrum of the electrons via a — 

(7e - l)/2. 

We define i?syn as the electron energy beyond which syn- 
chrotron and inverse Compton losses are significant, hence 
the ratio K(E) is not constant, and the equipartition estimate 
fails (see Sect. I3.4> . In the energy range Ep < E < E^y^, 
the ratio is constant (K = Kq) and the proton and elec- 
tron spectra are power laws with the same spectral index 
(7c = 7p = 7) if nonthermal bremsstrahlung loss of the elec- 
trons is negligible (see Sect. l3.4> . Then the radio synchrotron 
spectrum can be used to extrapolate the proton spectrum, so 
that Up in Eq. (IA13> can be replaced by: 



|(7-l)/2 



[c2h)E^lB 



7 7 r(7+1)/2- 



and the total cosmic ray energy density becomes: 



The projected field component B±^ has to be replaced by 
the total field B: 



B 



(7+l)/2 



(A17) 



C4(i) = [cos(i)]'^'''+^^/^ is valid for observation of a region 
where the field is completely regular and has a constant incli- 
nation i with respect to the sky plane (i = 0° is the face-on 
view), e.g. in a small region of the Milky Way or a galaxy. 
If the field is completely turbulent and has an isotropic an- 
gle distribution in three dimensions, B\ ~ {2/3)B^ and 
£4(1) = (2/3)^'''"'"^-'/^. If the synchrotron intensity is aver- 
aged over a large volume, [cos(i)]''''+^^/^ has to be replaced 
by its average over all occuring values of i. The case of a field 
oriented parallel to a galactic disk was treated by Segalovitz 
et al. (1976. their Eq. (3)). 

Combining ecR = cb = B'^^/Stt with Eq. ( IA16> and 
replacing 7 by (2a + 1) leads to: 



Beq = { 4^(2a + 1) (Ko + 1) L 4"'" (i^/2ci)" 



/ [(2a-l)c2(a)Zc4(z)]} 



l/(a+3) 



(A18) 



For use in the classical formula, Eq. ( IA17t has to be re- 
placed by: 



(A19) 



With the definitions introduced before, we can write the 
classical minimum-energy formula Q as: 



B, 



:iass = {67r(/C + l)/. {iy/2r q'/^ [i^^l'-'^^ 

-y'iLl / [(2« - 1) C2(a) I C4(*)] }'^'(A20) 

A table to compute -Bciass from ly (in mJy/arcsec^) and I (in 
kpc) is given in Govoni & Feretti ( 2004). 

As discussed in Sect.|2 the integration of the radio spec- 
trum between the fixed frequencies fmin and fmax introduces 
an implicit dependence on field strength B which is ignored 
when computing the derivative detot / dB. The constant expo- 
nent 2/7 in Eq. ( IA20t is the result of this error, which is one 
of the reasons why the classical formula gives wrong results 
(see Sect. 2). 



ecR = 7 (Ko + 1) E^--^ (v/2ci) 
/[2{^-2)c,{^)lB^2'''^''] 



(7-l)/2 



(A16) 



Note that, according to Eq. (IA6> . ecR includes protons up to 
the highest energies. 



